Effects of dipole-dipole interaction between cigar-like BECs of cold alkali atoms: 

Towards inverse square interactions 
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We study the dipole-diploe interaction between ultracold alkali Bose atoms which, in general, was 
neglected because of the strong short range coupling in the alkah atoms. The heavy suppression 
of the short range coupling and enhancement of the dipole-dipole coupling between the cigar-like 
Bose-Einstein condensates(BECs) imply that the latter becomes dominant. This gives a proposal 
to observe the effect of this dipole-dipole interaction. In the limit of long length of the single 
BEC, the resulting effective one-dimensional models possess an effective inverse square interacting 
potential, the Calogero-Sutherland potential, which plays a fundamental role in many fields of the 
contemporary physics but their direct experimental realization challenges the experimentalists for 
a long time. 

PACS numbers: 67.85.Pq, 37.10.Jk, ll.30.Pb 
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Introduction The dipole-dipole interactions between 
cold atoms/molecules play an important role in quan- 
tum simulations of strong correlated systems with cold 
atom/molecule gases lu |2|. Most of these studies are for 
the polar molecules [3|, the atoms with relevant strong 
dipolar couplings Q , and cold Rydberg atoms . The ef- 
fect of the dipole-dipole interaction between alkali atoms 
was thought of negligible because their dipolar coupling 
is too weak to compare with their s-wave scattering. 

Experimentally, the cold alkali atom gases are much 
easier to be made and manipulated, and are more sta- 
ble than those polar molecule and dipolar atom gases. 
An interesting question then is: Is there any way to 
change the strength of both kinds of interactions be- 
tween alkali atoms so that the dipole-dipole interaction 
becomes dominant or the s-wave scattering becomes neg- 
ligible? If so, we can use the alkali atom gases to study 
the dipole-diploe interactions, which have a great exper- 
imental advantage over using the other less stable cold 
atom/molecule systems. 

One way to reduce the s-wave scattering is tuning the 
s-wave interaction to zero using a magnetic Feshbach res- 
onance and then one can measure a very weak dipole- 
dipole interaction Q . Another way to reduce the s-wave 
scattering is to consider the effective interaction between 
the single mode fluctuations in an array of cigar-like cold 
alkali atom Bose-Einstein condensates (BECs)which are 
confined in a one-dimensional optical lattice @, 01 ■ How 
does the dipole-dipole interaction change for these single 
modes? This will be main topic in this paper. We will 
show an amazing result, i.e, the dipole-dipole interaction 
between atoms may reduce to an inverse square interac- 
tion between the single modes while the strength of the 
coupling constant is enhanced proportional to the atom 
number in a single BEC. On the other hand, the effective 
coupling of the s-wave scattering is reduced roughly as 
the inverse of the atom number in a single BEC. There- 
fore, in the limit of long length of the single BEC, the 





FIG. 1: (color online) The array of cigar-like cold atom clouds 
in an optical lattice, (a) The dipoles are polarized along 
the i-directions. (b) The dipoles are polarized along the z- 
directions. The yellow spots are the BEC substrate and the 
black are the quasi-atom with an inverse-square interaction. 



effective on-site interaction can be neglected and the sys- 
tem becomes a pure inverse square interacting one. 

One-dimensional inverse square interacting many-body 
systems, starting from Calogero-Sutherland (C-S) mod- 
els [1, nowadays play a fundamental role in many 
branches of modern physics [3 [HI 



Besides their in- 
terest and importance in theoretical and mathematical 
physics, the C-S-type models are often applied to ex- 
plain the physical phenomena in fractional quantum Hall 
effects 12 1, quantum exclusion statistics [13[, and black 
hole physics [3], and so on. However, the existence of 
the inverse square interaction lacks a direct experimen- 
tal realization. Our proposals in this paper present a 
possibility to realize them in laboratories. 

We consider a cold Bose atom cloud with dipolar in- 
teractions. With a very strong z-direction confinement, 
the cloud has a pancake-like shape trapped in the x-y 
plane with trapping frequencies cOx and ujy . Using a peri- 
odic optical potential along the x-direction, the pancake 
is incised into a lattice of cigar-like gases. An applied 
external field parallel to the x-axes, or perpendicular to 
this x-y plane polarizes all dipoles parallel to the external 
field. ( See Fig. [TJa) or (b), respectively). This setup is 
towards the Calogero modelQ. 

Another setup is rolling the x-y plane to the surface of 
a cylinder. This can be obtained by confining a magnetic 
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FIG. 2: (color online) The cigar-like BEG (green) are con- 
fined in a cylinder, (a) The magnetic dipoles; (b) The electric 
dipoles. The inserts indicate the directions of the dipoles. 



dipolar BEC cloud to the surface of the cylinder with 
a stable current passing through the axis or an electric 
dipolar BEC cloud with a charged bar placed in the axis. 
(See FigJ^a) or (b), respectively.) Sutherland's periodic 
variant [9| can be realized with this setup. A close rela- 
tive of such a setup, cold atoms in a ring-shaped optical 
lattice, has been experimentally realized [l^ . 

Reduce to a single-band model We first focus on the 
setup of Fig. [T]and set = 0. The energy scale in the 
problem is hujy ksT ^ hujp where Up is the frequency 
of lattice potential along the z-direction. Thus, we face 
a multi-band problem. Following the two-step procedure 
proposed in [7| , we reduce the multi-band problem to an 
effective single-band one. First, we treat a single cigar- 
like gas in a given lattice site. The lowest energy level is a 
BEC ^0(2/) and higher levels are the thermal cloud. The 
equation of state with only short range interactions has 
been exactly solved in the weakly interacting limit [61] . It 
shows that even if there are phase fluctuations (collective 
modes), the atom number in a BEC may be determined 
at every given temperature For the atom carrying a 
dipole, besides the short range interaction, two dipoles 
located in r and r' interact as 



d d' 3(d R)(d' R) 

i?3 



(1) 



where R = r — r' and R = R/i?; d and d' are the dipole 
moments of atoms located at r and r', respectively. We 
consider all d are parallel as shown in Fig. [TJ The Fourier 
components of the dipolar interacting potential do not 
contain any singularity in the momentum space but only 
provide an anisotropic term. Therefore, for a cold atom 
gas with weak dipolar interaction, we can also use the 
method given in Ref. [1] to solve the equation of state 
and determine the atom number No{T) in a BEC. 

The second step is to consider the coupling between 
sites If TT-ujy ^ NylJQ ^ hujp, where Ny is the av- 
erage number in a cigar-like gas and Uq the on-site re- 
pulsion, the wave function may be approximated by the 



product of the single atom ground state wave function 
in the a;-direction and the cigar-like atom gas in the y 
direction. In this approximation, the atomic field opera- 
tor may be written as ipir) = J2i v cLi,vW,^{x — Xi)'^^{y), 
where u is the band index; Wy(a;) is the Wannier function 
along the x-direction and ^,^{y) are the wave functions 
in the y-direction for the v-th band . a.i^^ is the corre- 
sponding annihilation operator for a Wannier state in the 
iz-th band at each site. 

Because we have a BEC ^'o(j/) at every site and iVo ^ 
Ny, the coupling between sites will be dominated by tun- 
nehng from a BEC to condensate as opposed to a BEC 
to thermal cloud. The multi-band problem may be re- 
duced to a single band one 0], i.e., the field operator 
may be approximated by ipir) « J^i'^i'^i^ ~ ^i)'^o{y), 
where the bosonic mode Oj = Uifi corresponds to an- 
nihilating a quasi-atom mode described by the product 
of the Wannier function w{x) = wq{x) and the cigar- 
like BEC wave function. With this approximation, we 
can reduce the multi-band problem to an effective single 
band boson model with both renormalized on-site and 
renormalized dipole-dipole interactions. The renormal- 
ized hopping tn and interaction Ur can be estimated 
in the Thomas- Fermi approximation [7]. The hopping t 
may not be renormalized [16j | while Ur is reduced from 
its bare value considerably due to the repulsive on-site 
interactions that spread out the condensate wave func- 
tion. The consistent condition for this reduction of Ur 
is Ip/a <^ Nq <^ {hujy/hujp)'^\/2TrLy/a, where a is the 
s-wave scattering length. Ip is the lattice spacing and 
hiijp = h^/mlp] and L^^y are defined by ^hjmujx.^y- 

The lattice model with an inverse square interaction 
The effective single-band lattice model can be obtained 
by keeping only the nearest neighbor hopping and ex- 
panding on-site energy to the second order near mean 
occupation Nq 

0,[I1, i.e.. 
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where Sui = alai — Nq is the deviation of the atom num- 
ber from the average number per site; t — J dx w* {x — 
Xi){—j^-^-\-Vp{x))w{x — Xj) for a pair of nearest neigh- 
bor sites. We emphasize that here the hopping 'particle' 
is not the whole cigar-like BEC as a 'giant particle' but 
only the particle number variance between the adjacent 
lattice sites. The renormalized on-site potential Ur, is 
given hy Uq^ ^ ^^^TT 1^1 • Hereafter, we take Ip ^ 1 
as the unit of length except explicitly stating. The renor- 
malized dipolar interaction potential is given by 



Ud, 



drdr'\w{x — Xi)\'^\w{x' — Xj)\^ 



Vd{r,r')\My)nMy')\ 
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We now assume the condensate size Ry (the Thomas- 
Fermi radius of the cigar-like BEC in the y-direction) is 
much larger than the lattice size L^. Then, one can use 
the approximation |^o(2/)P ~ the average density 
along in a single BEC and Ry ^ L^. Under these ap- 
proximations, it is easy to integrate over y and y' and 
arrives at 

Ud,ij~ dxdx'\w{x - Xi)\'^\w{x' - Xj)\'^ (4) 



GRy 



{x-x'yJ{x~x')^ + Rl i^^-^jr Rtj 

where G = -U^Nq for Fig. [^a) and G = 2(PNq for Fig. 
[ijb). Thus, we have a one-dimensional lattice model with 
the on-site and inverse square interactions between the 
single modes from the ground state. We see that the 
finite size effect of the cigar-like BEC can only supply a 
negligible Of-^) correction. 

Domination of the renormahzed dipole-dipole interac- 
tion The stability condition of a single BEC at a given 
site is given by edd = d'^/{3Uo) < 1 [17|. For the alkah 
metal atoms, d ~ I^b and Cdd ^ 1, e.g., Cdd = 0.007 for 
^''Rb. While the renormahzed on-site interaction Uu is 
reduced by a factor from Uq = dis enhanced 

to djf = d\/No as shown in (jH). Thus, although edd ^ 1, 
the renormahzed on-site interaction may be much weaker 
than the renormahzed dipole-diploe interaction, i.e., the 
renormahzed ed^da = <Pr/{Wr) = eddNoRy/lp-, may be 
much larger than the unit. For example, if A^o ^ 10'^ and 
Ry/lp > 50 for ®^Rb, ed^dR Z 10^. Therefore, in the ef- 
fective one-dimensional single-band model, the renormal- 
ized dipole-dipole interaction dominates and the renor- 
mahzed on-site interaction may be neglected. 

Trapping potential in the x-direction We have obtained 
a homogeneous lattice boson model with inverse square 
interaction by neglecting the trapping in the x-direction. 
Adding back the harmonic trap with the frequency uj^ 
may affect the atom number Nq{x) at each site. If 
the width of a single BEC is wider, the number A'o of 
atom per site may vary site by site and Aq in the cen- 
ter of trapping may be several times of Aq in far away 
from the center. Within a cigar-like single dipolar BEC, 
the dipole-dipole repulsion is strong enough against rais- 
ing the width of the single BEC at the trapping center. 
Therefore, the factorization of the ground state wave is 
still valid and Aq keeps an constant. The resulting inho- 
mogeneous model is _ff in ([2]) added by a term Vh.i^ni 
where 



1 



Vh,i = T^'^i^x I dx\w(x - Xi)\ 



2 2 

X 



(5) 



is the harmonic trapping potential along the a;-direction. 

Dilute gas limit and C-S model Assuming the fluctua- 
tion atom number 5N = 5ni is fixed. In the dilute gas 



limit with 5N/Lx <C 1, the JA'-particle system may be 
described by a continuous model because the dispersion 
—tcosk/h ^ A;^/(2m/j) with to^ ^ h^/t. Neglecting the 
on-site interaction as argued above, the effective Hamil- 
tonian reads 



SN 



2mR dxf 



ruRUj^ ^ 



G 



(6) 



where 



Xk and the renormahzed trapping po- 



tential is defined by tow^ = mRLo\. The model described 
by the Hamiltonian © is the famous C-S model [1, . 
A bosonic 'particle' at Xi in this continuous model cor- 
responds to one more atom than Aq at the lattice site i 
in the lattice model. The ground state wave function is 
given by 



^q^x{xi,...,xsn) 



n 

l<3<k<SN 



Xjk\ G 



2h / , j 



(J) 



with the ground state energy Eg — ^5NTiujr{\{5N — 
1) -I- 1). Here, A is the solution of the equation A(A — 1) = 
{mR/n^)G, e.g., g ^miNoUo/t for the repulsive 
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interacting ^'^Rb atoms. This gives A± = i±v^±M, por 
a weak on-site interaction with U^/t < 1 ov \, g may 
vary as t and is of the order 0(1) and so is A if A'o is 
fixed to the order of 10"^. 

For the setups in Fig. 1(a), g < 0. If 5 < -1/4, 
A becomes imaginary. This gives an imaginary ground 
state energy and the system is not stable against the 
collapse. However, when — 1/4 < 5 < 0, the ground state 
is stable with < A_ < 1/2 although the interaction 
term is formally negative. 

For the setups in Fig. 1(b), g > 0. For < A_ < 0, 
i.e, 5 < |, the wave function is still square integrable. 
This implies if the repulsive between particles is not 
strong enough, the system will also collapse. For a re- 
pulsion with g > I , ([7]) is no longer square integrable 
for A = A- < — ^ and a physically stable solution is given 
by A = A+ which is larger than 1 here. 

Periodic interaction potential We now consider the ge- 
ometry in Fig. [5J In this geometry, the cylinder extends 
along the ^/-direction and circles in the x-direction. If 
Ry 3> one may prove that, up to a constant terms, 
the renormahzed dipole-dipole interaction is a periodic 
inverse square interaction 



7r^A(A - 1) 
Li sm^l-KXij / Lx] 



(8) 



where a_ = — 2 and a+ = 1 with respect to Fig. 

and for Fig. [Ifb). The similar but periodic solutions 

were obtained by Sutherland |9|. 

Experimental implications The first effect from the 
dipole-dipole interaction between the alkali atoms is the 
collapse of the single modes to a single lattice site (a given 
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FIG. 3: (color online) (a) The schematic diagram of the mo- 
mentum distribution, (b) The top-view of the time-of-flight 
image. The grey area is the original cylinder. The yellow only 
has very few atoms. Most of atoms arrive at the time t in the 
red area and the green is of those atoms with the momentum 
distribution (fTU)) . 



single cigar- like BEC), when g < —1/4 for the setups of 
Fig. [l]^a) and Fig. ^a) and g < 3/4 for the setups of 
Fig. [TJb) and Fig. [UJb). This will never happen for 
a system dominated by the on-site repulsive interaction. 
The time-of-flight experiments measures the momentum 
distribution of the trapped cloud 



18|. In these experi- 



ments, this collapse can be observed and can be thought 
as a definitive evidence of the effect of the renormalized 
dipole-dipole interaction. 

In addition to this definitive evidence for the effect 
of the dipole-dipole interaction in alkali atoms, one can 
measure more observable of the systems. For example, 
one can consider the long wave length limit of the Green's 
function at t = of the Sutherland model, which is given 
by [19| 



easy to be distinguish from the cae with no dipolar in- 
teraction. However, for A > 2, we switch off all traps but 
the inner wall of the cylinder. Because the cloud along 
the radial direction is the most squeezed, its expansion 
along this direction is much faster than that along the 
axial and tangent directions, owing to a larger zero point 
energy. Assume uq is the velocity of the radial motion 
of a single cigar-like cloud and Vm is the fastest velocity 
of the particles along the circle. Let us see the top view 
of the cylinder. At the time t after turning off the traps, 
most of background BEC atoms arrive at the red area in 
Fig. IS^b) while the momentum distribution in Eq. (fTO)) 
can be figured out from the atoms arrive at the green . 

To match the theoretical results with the experimental 
data more precisely, the Green's function ([9|) and the 
momentum distribution (llOp should be calculated in term 
of the lattice model ([2]). However, we believe Eq. ([TU| is 
qualitatively correct if G is in the designed region. 

In conclusions, we showed the weak dipole-dipole inter- 
action between the alkali atoms may renormalize to an 
inverse square potential dominating over the renormal- 
ized on-site interaction even without changing the bare 
on-site interaction by any external field. On the other 
hand, our proposal raised the possibility to implement 
the C-S model in nature, which was long expected but 
not directly realized yet in laboratories. 

The author thanks Yingmei Liu, Kun Yang, Su Yi and 
Li You for useful discussions. This work was supported 
by National Natural Science Foundation of China, the 
national program for basic research of MOST of China. 
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G(x,0) cx l + 2^{-iy"B. 



cos 2mkpx 



(9) 



where Bm are regularization-dependent constants and 
kp — -kSN/ Lx is the effective Fermi momentum. The 
momentum distribution is the Fourier transformation of 
this Green's function. For the lowest energy sector with 
m — 1, the momentum distribution for small k and large 
k is approximated by (loj 



rifc oc |fc|^/2"^ if |fc| < fcF or cx |fc|" 



2A-2 



if |fc| > fci.(10) 



We see that when A < 2, the system is in a quasi- 
condensed state while the momentum distribution van- 
ishes in a power law for small k when A > 2, a Luttinger 
liquid behavior. At the critical point A = 2, a loga- 
rithmic divergence exists, i.e., Uk = {l/2)ln{2kp/\k\) for 
|fc| < 2kp and nj. = for |fc| > 2kp Graphically, 
these momentum distributions have been shown in Ref. 
[igj and we give them schematically in Fig. [Ha). 

The projected density profile in the plane perpendic- 
ular to the y-axis may directly reflect the momentum 
distribution Uk since the momentum along the circle is 
a good quantum number. The quasi-condensate is not 
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